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We review the notion of symmetry breaking and restoration within the frame of nuclear 
energy density functional methods. We focus on key differences between wave-function- 
and energy-functional-based methods. In particular, we point to difficulties encountered 
within the energy functional framework and discuss new potential constraints on the un- 
derlying energy density functional that could make the restoration of broken symmetries 
better formulated within such a formalism. We refer to Ref. d for details. 

1. Introduction 

Symmetries are essential features of quantal systems as they characterize their en- 
ergetics and provide transition matrix elements of operators with specific selection 
rules. However, certain emergent phenomena relate to the spontaneous breaking of 
those symmetries In nuclear systems, such spontaneously-broken symmetries (i) 
relate to specific features of the inter-particle interactions, (ii) characterize internal 
correlations and (iii) leave clear fingerprints in the excitation spectrum of the sys- 
tem. In finite systems though, quantum fluctuations cannot be ignored such that 
the concept of spontaneous symmetry breaking is only an intermediate description 
that arises within certain approximations. Eventually, symmetries must be restored 
to achieve a complete description of the system. 

2. Wave-function-based methods vs EDF-based method 

In wave-function-based methods, the symmetry breaking step, e.g. the symmetry 
unrestricted Hartree-Fock-Bogoliubov approximation, relies on minimizing the av- 
erage value of the Hamiltonian for a trial wave-function that does not carry good 
quantum numbers, i.e. which mixes irreducible representations of the symmetry 
group of interest. Restoring symmetries amounts to using an enriched trial wave- 
function that does carry good quantum numbers. One typical approach is to project 
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out from the symmetry-breaking trial state the component that belongs to the 
intended irreducible representation. Wave-function-based projection methods and 
their variants are well formulated quantum mechanically^. On the other hand it is 
of importance to analyze further their Energy Density Functional (EDF) counter- 
parts which have been empirically adapted from the former to deal quantitatively 
with properties of nuclei 

The single- reference (SR) EDF method relies on computing the analog to the 
symmetry-unrestricted Hartree-Fock-Bogoliubov energy £ as an a priori general 
functional £[p" , k" , n qq *] of the density matrices p" and k" computed from the 
product state \<&(g)). Here, the label g denotes the parameter(s) of the symmetry 
group Q of interest, e.g. Euler angles for 5*0(3) or the gauge angle for f/(l). The 
multi-reference (MR) EDF method, that amounts to restore symmetries broken 
at the SR level, relies on computing the analog to the non-diagonal energy kernel 
£[g',g] associated with a pair of product states (®{g')\ and | $(<?)) as a general 
functional of the transition one-body matrices, i.e. £[g',g] = £[p 9 9 , k 9 9 , k" *]. 
From such a energy kernel, the symmetry-restored energies £ x are obtained through 
an expansion ^ over the volume of the symmetry groupEI 



£[g\g]N[g' ,g]^Y, c ^ c ^ £X S ab(9-9') , (1) 

Xab 

N[g\g}^J2 c ^ox b S x b (g~g / ) , (2) 



Xab 



where N[g',g] = ($((/)|$(g)), while S^ b (g—g') denotes a unitary irreducible repre- 
sentation (labeled by A) of the symmetry group. 

A key point is that, as opposed to wave-function-based approaches, the EDF 
method does not rely on computing £[g' ,g\ from the average value of a genuine 
scalar operator H such that all expected mathematical properties of such a kernel 
are not necessarily ensured a priori. Consequently, one may wonder whether the 
symmetry constraints imposed on the energy kernel £ [p, k, k*] at the SR level SI are 
sufficient or not to making the MR-EDF method well defined? As a matter of fact, 
a specific set of constraints to be imposed on the non-diagonal kernel £[g' , g] have 
been work ed out t o fulfill basic properties and internal consistency requirements^. 
Still, Refs. I 6 | 7 | 8 | 9 | fr ave s hown in the case of the U(l) group, i.e. for particle-number 
restoration (PNR), that such constraints were n ot su fficient to making the theory 
well behaved. In particular, it was demonstrated that Fourier components £ N 
could be different from zero for a negative number of particles N < 0. Contrarily, it 
can be shown that E N is zero ^ for N < when it is computed as the average value 
of a genuine operator H in a projected wave-function, i.e. in the wave-function- 
based method. Applying the regularization method proposed in Ref. the nullity 
of the non-physical Fourier components was recovered^ 



a Such an expansion is nothing but the Fourier decomposition in the case of the U(l) group 
associated with particle-number conservation. 
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3. Towards new constraints? 

The case of U{1) was particularly instructive given that clear-cut physical argu- 
ments could be used to assess that certain coefficients of the (Fourier) expansion 
of the energy kernel should be strictly zero. Such an investigation demonstrated 
that the MR-EDF method, as performed so far, faces the danger to be ill-defined 
and that new constraints on the energy kernel £[g',g] must be worked out in or- 
der to make the symmetry-restoration method physically sound. The regularization 
method proposed in Ref. ^ that restores the validity of PNR can only be applied if 
the EDF kernel £ [p, k, k*] depends strictly on integer powers of the density matri- 
ces^, which is an example of such a new constraint. 

For an arbitrary symmetry group, the situation might not be as transparent as 
for U(l). Indeed, it is unlikely in general that certain coefficients of the expansion of 
£[<?', g]N[g' , g] over irreducible representations of the group are zero based on phys- 
ical arguments. The challenge we face can be formulated in the following way: al- 
though expansion [1] underlining the MR-EDF method is sound from a group-theory 
point of view, additional mathematical properties deduced from a wave-function- 
based method must be worked out and imposed on £ [g',g] to make the expansion 
coefficient £ x physically sound. The next section briefly discusses an example of 
such a property in the case of SO(3), i.e. for angular momentum restoration, that 
could be used to constrain the form of the functional kernel £[fi',Q]. More details 
can be found in Ref. H. 



3.1. Mathematical property associated with angular-momentum 
conservation 

We omit spin and isospin for simplicity and consider the rotationally-invariant 
nuclear Hamiltonian H = T + V in which three-nucleon and higher many-body 
forces are disregarded for simplicity. Considering an eigenstate |0 LM ) of I? and L z , 
as well as using center of mass R = (r*i + r2)/2 and relative coordinates r = r± — ?2, 
tedious but straightforward calculations allow one to show that the potential energy 
takes the form^ 

^ S % M LM^L L M ) ) = jdRV^{R) , (3) 

where the potential energy density V LM (R) thus defines from the two-body density 

, . [2] LMLM i 

matrix p - read 
1 Rf 



V LM {R) = £ Cxmi/o vf} L , (R) Y°, (R) . (4) 



L'=0 



The weight vf} L ,(R) depends on the norm of R only and is related to a reduced 
matrix clement of the two-body density matrix operator recoupled to a total angular 
momentum L' . The remarkable mathematical property identified through Eq. |4] 
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is that the scalar potential energy V L is obtained from an intermediate energy 
density V LM (R) whose dependence on the orientation of R is tightly constrained 
by the angular-momentum quantum number of the underlying many-body state 
|O iM ), i.e. its expansion over spherical harmonics is limited to L' < 2L. Such 
a result is unchanged when adding the kinetic energy (density) to the potential 
energy (density) such that we restrict ourselves to the latter for simplicity. Of 
course, the energy eventually extracts the coefficient of the lowest harmonic, i.e. 
V L = JdRv [ H(R). 

3.2. Wave-function-based versus EDF-based methods 

Since property [4] is general, it is straightforward to show that it can be recovered 
within the frame of the wave-function-based symmetry- restoration method 

The key point to underline relates to the fact that property [4] cannot be de- 
rived a priori in the EDF method. Indeed, the potential energy part of the kernel 
£[fi',fi] is not explicitly related to the two-body density matrix in such a case. 
Consequently, there is no reason a priori that the energy density V LM (R) displays 
property HI i.e. the angular dependence of V LM (R) is likely to display harmonics 
Y®,(R) with V > 2L. One might argue that it is not an issue considering that the 
symmetry-restored potential energy V L eventually relates to the harmonic Yq(R) 
only. However, a formalism that provides V LM (R) with a spurious angular con- 
tent will certainly also provide the coefficient Vlo(R) of the lowest harmonic with 
unphysical contributions. To state it differently, it is likely that constraining the 
MR-EDF kernels £ [p n ' n , n n ' n , n nn ' *] to produce an energy density £ LM (R) that 
fulfils the mathematical property |4] will impact at the same time the value of the 
weight £lo(R), and thus the value of £ L . To some extent, this is similar to the 
situation encountered with U(l) where restoring the mathematical property that 
Fourier coefficients £ N with N < should be strictly zero impacted the value of all 
other Fourier coefficients^. 



4. Conclusions 

We briefly review the notion of symmetry breaking and restoration within the frame 
of nuclear energy density functional (EDF) methods. Multi-reference (MR) EDF 
calculations are nowadays routinely applied with the aim of including long-range 
correlations associated with large- amplitude collective motions that are difficult 
to incorporate in a more traditional single- reference (SR), i.e. " mean- field" , EDF 
formalism ^^1. 

The framework for MR-EDF calculations was originally set-up by analogy with 
projection techniques and the Generator Coordinate Method (GCM), which are 
rigorously formulated only within a Hamiltonian/ wave- function-based formalism^. 
We presently elaborate on key differences between wave- function- and energy- 
functional-based methods. In particular, we point to difficulties encountered to for- 
mulate symmetry restoration within the energy functional approach. The analysis 
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performed in Ref.^to tackle problems encountered in Refs.D3II2IS]f or particle num- 
ber restoration serves as a baseline. Reaching out to angular-momentum restora- 
tion, we identify in a wave-function-based framework a mathematical property of 
the energy density E LM (R) associated with angular momentum conservation that 
could be used to constrain EDF kernels. Consequently, possible future routes to 
better formulate symmetry restorations in EDF-based methods could encompass 
the following points. 

(1) The fingerprints left on the energy density E LM (R) by angular momentum 

conservation in a wave-function-based method could be exploited to constrain 

the functional form of the basic EDF kernel £\p, k, k*]. 

171 

(2) The regularization method proposed in Ref.^to deal with specific spurious fea- 
tures of MR-EDF calculations should be investigated as to what impact it has 
on properties of the energy density £ LM (R) in the case of angular momentum 
restoration. 

(3) Similar mathematical properties extracted from a wave-function-based method 
could be worked out for other symmetry groups of interest and used to constrain 
EDF kernels. 

Efforts in those directions are currently being made. 
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